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A fixed-bias spherical Langmuir probe is included as part of the Vector Electric Field Instru- 
ment (VEFI) suite on the Communication/Navigation Outage Forecast System (C/NOFS) satellite. 

C/NOFS gathers data in the equatorial ionosphere between 400 and 860 km, where the primary 
constituent ions are H + and 0 + . The ion current collected by the probe surface per unit plasma 
density is found to be a strong function of ion composition. The calibration of the collected cur- 
rent to an absolute density is discussed, and the performance of the spherical probe is compared to 
other in situ instruments on board the C/NOFS satellite. The application of the calibration is dis- 
cussed with respect to future fixed-bias probes; in particular, it is demonstrated that some density 
fluctuations will be suppressed in the collected current if the plasma composition rapidly changes 


along with density. This is illustrated in the observation 
[http://dx.doi.org/10.1063/L4766333] 

I. INTRODUCTION 

A commonly used tool for plasma diagnostics in the 
upper atmosphere is the Langmuir probe. 1 The technique 
typically consists of sweeping the bias of a small metal- 
lic probe (typically a sphere or a cylinder) relative to the 
plasma potential to infer such parameters as plasma density 
and temperature. 2 Langmuir probes have flown on numerous 
satellite missions for the diagnosis of space plasmas, includ- 
ing Atmospheric Explorer,’ Dynamics Explorer 2, 4 and the 
Pioneer Venus Orbiter. 5 Another common application is to 
maintain a fixed bias on the probe in either the ion or electron 
saturation regime in order to provide a rapid measurement of 
density fluctuations. 6 This paper will focus on the ion satura- 
tion mode, where a negative bias is applied to ensure that only 
ions are collected. 

The Langmuir probe on the Communication/Navigation 
Outage Forecast System (C/NOFS) satellite, described in 
Sec. II, is a fixed-bias probe that measures the total ion den- 
sity in the topside equatorial ionosphere. This region typically 
consists of a mixture of 0 + and H + ions (and to a lesser ex- 
tent, He + and N + ). The transition height between 0 + and 
lighter ions has been studied extensively as a function of so- 
lar activity. 7 The nightside transition height near the equator 
is found around 600 km during a typical solar minimum, but 
can be lower than 500 km during the extreme solar minimum 
of 2008. 8 During a typical solar maximum (FI 0.7 = 200 sfu), 
this height can be well above 800 km. 9 Additionally, the day- 
to-day variation of ionospheric composition in the range of 
the C/NOFS satellite is significant. 1 ! The effect of ion com- 
position on the expected collected current is discussed in 
Sec. III. 

A calibration technique for a small spherical probe in a 
cold collisionless plasma where composition is known is pre- 
sented in Sec. IV, along with the application of this technique 
to the probe on C/NOFS. The performance of the calibration 

^Electronic mail: jeffrey.klenzing@nasa.gov. 

0034-6748/2012/83(1 1)/1 14501 /5/$30. 00 83, IF 


of plasma density enhancements on C/NOFS. 


is discussed in Sec. V, as well as the implication for future 
probe designs. 

II. INSTRUMENT DESCRIPTION 

The C/NOFS satellite was launched in April 2008 into 
a 13° inclination orbit with perigee near 400 km and apogee 
near 860 km. C/NOFS is equipped with multiple instrument 
suites designed to study the ion and neutral populations and 
their effect on the propagation of communication signals. 11 
The main focus of this paper will be the fixed-bias Langmuir 
probe (referred to as the Trigger Probe, or TP), which is part 
of the Vector Electric Field Instrument (VEFI) suite. 12 

The trigger probe on the VEFI is a fixed-bias partial- 
spherical Langmuir probe designed to measure relative 
plasma density fluctuations, as seen in Figure 1. The probe 
consists of a collector and a guard element which together 
form a sphere with a radius of 1 .25 cm. The guard is included 
to shield the bulk of the collector from effects of the support 
boom breaking the spherical symmetry and is held at the same 
potential as the collector. The collector is a partial sphere 
with a half-angle of 111.2° and has a total surface area of 
13.35 cm 2 . Both the collector and the guard are coated in TiN. 
The probe is mounted so that the support boom is perpendic- 
ular to the satellite ram velocity. The boom is 0.5 m long. The 
probe can be set to one of the four bias settings: — 1 .25, —2.5, 
—3.75, and —5 V. For these settings, the probe is well within 
the ion saturation region. The normal operating mode of the 
probe is —3.75 V. 

The probe is primarily used as a trigger for the burst 
memory on board C/NOFS to save high-resolution (>2000 
samples/s) electric field and density data. Nominally the probe 
collects data at 16 samples/s. As the C/NOFS satellite has 
been operational for more than four years, a large statistical 
database of ion densities is available for study. 

Several other instruments on board C/NOFS are em- 
ployed in the calibration procedure. The Planar Langmuir 
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FIG. 1 . The trigger probe (TP) for the VEFI suite on the C/NOFS satellite. 
The collector and guard form a spherical probe, with the guard removing the 
non-uniformities introduced by the boom. 


Probe (PLP) built by the Air Force Research Laboratory 
(AFRL) consists of a swept Langmuir probe and an ion trap. 1 3 
The ion trap is used to provide absolute densities for calibra- 
tion. Additionally, the Ion Velocity Meter (IVM), built by the 
University of Texas at Dallas, provides relative ion compo- 
sition data using the well-known retarding potential analysis 
technique. 14 


III. PROBE THEORY 

The equations describing the current collected by vari- 
ous probe geometries for different ion distributions have been 
discussed at length elsewhere. 13 ' 16 Here we will follow the 
technique of Hoegy and Brace. 16 The ion current collected by 
a sphere immersed in a thermal (Maxwellian) plasma consist- 
ing of a single species of ion is 


/ random — A probe 77 1 


I k B Tj 
2tt m j ’ 


(1) 


where qo is the fundamental charge, A prohe is the surface area 
of the sphere, Nj is the total ion density, k B is the Boltzmann 
constant, 7) is the ion temperature, and in, is the ion mass. If 
the sphere is biased relative to the plasma (at (S > app ) or moving 
at velocity uq relative to the plasma, the modified current can 


be described by 


/ sphere — / random 
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r 2 


erf(r) 


where 


up 

y/2k B Ti/mi ’ 


( 2 ) 

(3) 


= 


qo^fipp 

k B Tj 
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Note that as the velocity (uq) approaches zero, the error 
function can be approximated by a Taylor expansion, such 
that 

2 

■s/n 

(5) 

Equation (2) then reduces to the usual orbital motion limited 
collection equation 

/ sphere ~ 7 random ( 1 T 1] ) • (6) 

For the C/NOFS satellite, the ion temperatures encountered 
during the extreme solar minimum are between 600 K and 
1300 K 17 and slowly increasing through 2011 as solar activity 
rises. The spacecraft velocity is ~7.5 km/s, thus expected val- 
ues of r are ~l-2 for H + and >6 for 0 + . The probe is biased 
negatively so that the electron current is negligible. Further 
discussion of the ion collection for other cases (such as low 
Mach numbers and probes near the floating potential) can be 
found in the works of Hutchinson. 18-20 

For a plasma consisting of H + and 0 + , the average cur- 
rent to the sphere per unit density can be represented as 

Sphere = = /h+ A t+ + fo+Io+, (7) 

where /h+ and Io+ represent the average current contributions 
per unit density of H + and 0 + , and fo+ is the fraction of 0 + . 
Figure 2 shows these average current contributions as calcu- 
lated by Eq. (2) as a function of ion temperature and refer- 
ence potential for typically expected values in the range of the 
C/NOFS satellite. While there is a small effect in the expected 
current per unit density for fluctuations in ion temperature (a) 


erfir) 1 2 
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FIG. 2. The variation of expected current collected per unit density as a function of (a) ion temperature, (b) reference potential, and (c) ion composition. 
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or spacecraft potential (b), these changes are commensurate 
with an uncertainty of a few percent in the total 0 + composi- 
tion (c). This is because the currents due to H + and 0 + differ 
by more than a factor of six. Note that for the stationary unbi- 
ased sphere of Eq. (1), the difference in currents is expected 
to be only a factor of four (the square root of the ratios of 
masses). The additional current from H + is a result of both 
the negative bias of the probe and the relative speed between 
the probe and the plasma. 

For a typical ionospheric mission, 0 + will be collected 
only on the ram surface of the probe, while light ions such 
as H + and He + are collected over the whole surface. This is 
due in part to the large thermal velocity of light ions, but it is 
mostly due to the effect of the negatively biased probe on the 
motions of the ions. Light ions that are not in the direct path 
of the probe are more easily accelerated toward the probe and 
collected on all surfaces, as demonstrated by the simulations 
of Seran et al . 21 Heavier ions such as 0 + are accelerated to- 
ward the sphere as well, but the trajectories are not changed 
enough for these ions to be collected on the anti-ram surface. 
This is illustrated in Figures 8 and 9 of Seran et al. Addi- 
tionally, it is possible under the right combination of condi- 
tions that light ions may be collected in greater abundance on 
the anti-ram surface of a spherical probe than on the ram sur- 
face due to this acceleration effect (see Figure 1 1 of the same 
paper). 


TABLE I. Data used for each fit. 


LPBias (V) 

Days 

Number of points 

-1.25 

105-174, 2009 

79 295 

-2.50 

230-299, 2009 

138 810 

-3.75 

350, 2009; 54, 2010 

140 006 

-5.00 

155-179, 2008; 55-84, 2010 

83 340 


The calibration factor / is then fit to the function in 
Eq. (10) using a Levenburg-Marquardt 22 least-squares fit, 


X(/o+) = 


1 

a + bfo+ 


(10) 


The results of the fitting procedure for each bias setting (re- 
ferred to as LPBias) can be seen in Figure 3, where all data 
are plotted as a function of fo+. The mean and standard de- 
viations are plotted as error bars over the data points, and the 
fits are plotted as lines. The fit parameters a and b for all four 
bias settings are listed in Table II along with the correspond- 
ing coefficient of determination (R 2 ). All fits are performed 
using the MATLAB R2012a curve fitting toolbox. 

Figure 4 shows the mean and standard deviations for / 
from Figure 3 replotted as a function of applied potential 
(LPBias) for selected compositions. The changes in the mea- 
sured value of x for each bias setting are as expected. As the 


IV. CALIBRATION TECHNIQUE 

Because most of the expected variations in I sp here come 
from the plasma composition for the range of temperatures 
encountered by C/NOFS, we can rewrite Eq. (7) as a linear 
function of the fraction of 0 + : 


Isphere{fo + ) = ^H+ + (7o+ — 7h+)/o+- (8) 

The calibration of the VEFI trigger probe utilizes a multiplica- 
tive conversion factor x , defined as 

V i — Z trigger > (9) 

where logger is the current collected in nA and N; is the plasma 
density per cubic centimeter. Note that the conversion factor 
X is essentially the inverse of the effective collected current 
per unit density ( I sp here )■ The empirical fits described below 
are calibrated to the total plasma density provided by the PLP 
ion trap and the composition information as provided by the 
IVM. 

The value of / is determined by taking the ratio of the TP 
current to the total density as provided by PLP. For three of 
the four bias settings, 70 consecutive days of data are chosen 
to give ample coverage over all altitudes and local times (the 
precession of the C/NOFS perigee with respect to local time 
is about 65 days). The lowest bias setting (—5 V) only uses 
26 consecutive days in 2008 and 29 consecutive days in 2010 
due to the minimal time spent at this setting. The data are 
decimated so that one point is used every 64 s. The number of 
points for each fit is reported in Table I. 
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FIG. 3. The ratio of PLP density (V,) to trigger probe current (7tp) as a 
function of the fraction of O" 1 " as measured by IVM for the four bias settings. 
The cyan lines represent the mean and standard deviation of the data; the red 
lines are the fit to the data. 
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TABLE II. Fitted parameters. 


LPBias (V) 

a 

b 

R 2 

-1.25 

0.002978 

-0.002371 

0.9478 

-2.50 

0.004055 

-0.00331 

0.964 

-3.75 

0.005774 

- 0.004866 

0.9638 

-5.00 

0.007402 

- 0.006461 

0.9748 


applied bias becomes more negative, more ions are attracted 
to the probe and the resultant current per unit density in- 
creases. This leads to an overall reduction in x- The spread 
in x also decreases as the applied bias becomes large relative 
to any fluctuations in spacecraft potential. 


V. PERFORMANCE IN EARTH’S IONOSPHERE 

The performance of the calibrated fits is compared for all 
data in 201 1 in Figure 5. In addition to the calibration method 
described above, two other “densities” are calculated using 
an assumption that the probe is always in an environment con- 
sisting of a single ion species (H + or 0 + ). The vertical dashed 
lines represent a discrepancy of 30% between VEFI and PLP. 
The calibration utilizing the exact composition information 
typically provides densities that agree within 30% of the PLP 
density measurement, while the other two estimates of density 
may be off by as much as a factor of ten. 

The implication of this work is not just that ion com- 
position is vital to calculating an absolute density from a 
fixed-bias spherical probe, but also for calculating relative 
density fluctuations for satellites near the H + /0 + transition 
height. To illustrate this, Figure 6 shows a single orbit where 
the C/NOFS satellite passes through a series of plasma den- 
sity enhancements. These structures are frequently encoun- 
tered by C/NOFS in the nighttime ionosphere and consist of 
a sharp increase in both the total density and the relative frac- 
tion of 0+. 2 ’ Panel (a) shows the current measured by the 
trigger probe during a typical encounter with these structures 



LPBias (V) 


FIG. 4. The ratio of PLP density (Nf) to trigger probe current (/tp) as a 
function of the fraction of applied potential for four different levels of 0 + 
composition. 



FIG. 5. Performance of the calibration method for total density for all days 
in 201 1. The vertical dashed lines represent discrepancies of 30% versus the 
PLP density. 


on the night of September 15, 2009; panel (c) shows the rel- 
ative composition of 0 + . In panel (b), the absolute densities 
reported by both PLP and the calibrated TP are reported as 
solid lines, with a third option using a constant assumption of 
70% 0+ to convert the TP currents to densities (referred to as 
the uncalibrated VEFI densities). The uncalibrated curve sup- 
presses the density increase, leading to a loss of information. 
This is due to the increase in 0 + inside the plasma enhance- 
ment, which is collected by a spherical probe less efficiently 
that the H + outside the enhancement. In order to capture these 
density fluctuations, knowledge of the plasma composition is 
needed in addition to the collected current. 

This suppression effect is the result of flying an ion sat- 
uration probe through a two-species plasma where the two 


Orbit 7675 



FIG. 6. An example of plasma density structures where a changing compo- 
sition results in a suppression of the features if the probe is uncalibrated. 
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ion species have very different Mach numbers (r). For a high- 
velocity satellite (in this case, uq — 7.5 km/s), this results 
in Mach numbers of roughly 1 and 6 for H + and 0 + , re- 
spectively. Because of the large spacecraft velocity, varia- 
tions due to ion temperature are expected to be minimal (see 
Figure 2(a)). For a lower-velocity platform such as a sound- 
ing rocket (where spacecraft velocity is typically less than 
1 km/s), this is not the case. For sounding rockets, the flow 
of H + and 0 + will typically be subsonic (r < 1), and now 
fluctuations in ion temperature are expected to have a signifi- 
cant effect on the effective current collected. A suppression of 
density fluctuations similar to that observed on C/NOFS may 
occur for heavier species such as Fe + , where r may be 2-3. 

VI. SUMMARY 

The calibration technique for the spherical fixed-bias 
Langmuir probe for the VEFI on C/NOFS is described. It is 
found that the highly-variable ion composition in the vicin- 
ity of the satellite strongly affects the current collected by 
a spherical probe for a given density. The calibration of the 
probe is greatly improved by incorporating in situ measure- 
ments of ion composition as provided by the IVM instru- 
ment. Additionally, there are events such as plasma density 
enhancements where the changes in total density and relative 
ion composition combine in such a way that the probe cur- 
rent is relatively unchanged while the ion density doubles in 
magnitude. This is an important consideration for similar ion 
saturation probes designed to look at relative fluctuations near 
the H + /0 + transition height. For probes well above this tran- 
sition height (and therefore in a single species plasma), rela- 
tive density fluctuations can be reliably determined by a fixed- 
bias probe. However, knowledge of the in situ composition 
is essential in truly calibrating a probe in a multi-component 
plasma. 
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